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Abstract 

Using a density matrix description in space we study the evolution of wavepackets in a fluctuating 
space-time background. We assume that space-time fluctuations manifest as classical fluctuations of 
the metric. From the non-relativistic limit of a non-minimally coupled Klein-Gordon equation we 
derive a Schrodinger equation with an additive gaussian random potential. This is transformed into 
an effective master equation for the density matrix. The solutions of this master equation allow to 
study the dynamics of wavepackets in a fluctuating space-time, depending on the fluctuation scenario. 
We show how different scenarios alter the diffusion properties of wavepackets. 

1 Introduction 

The search for a quantum theory of gravity is still a work in progress and leaves many open questions. 
Currently there is no final version of a quantum gravity theory but it is expected that one of its con- 
sequences is the appearance of some kind of space-time foam. This can be thought of as space-time 
fluctuations which could manifest classically in the low-energy limit as fluctuations of the metric or of the 
connection. 

In the most simple version one can think of a Minkowskian background on which small space-time 
dependent metrical fluctuations are imposed. In this context it has been shown for quantum mechanical 
systems that such a scenario has the following consequences: (i) we obtain a modified inertial mass and, 
thus, an apparent violation of the weak equivalence principle |T1[51[S]. (ii) quantum systems will suffer 
decoherence induced by such quantum gravity induced space-time fiuctuations which were discussed in 
[4]- [7]. In the context of the propagation of light (iii) fluctuating light cones and (iv) angular blurring were 
discussed in [5] . Since space-time fluctuations can also lead effectively to non-localities they can manifest 
themselves as (v) modified dispersion relations [5]-[T^. 

(vi) higher-order time derivatives can appear in the context of equations of motions of higher order 
[T5] leading to non-localities w.r.t. the time variable. It might be expected that these terms appear as a 
consequence of quantum gravity scenarios. 

In this paper we investigate a further consequence of the model worked out in [3j. We consider the 
effect of space-time fluctuations on the evolution of wavepackets and calculate their modified evolution 
in terms of the mean-squared displacement and higher-order moments. The modification depends on 
the fluctuation scenario characterized by a spatial correlation function. Such kind of evolution equations 
have been discussed in |14j-[16j in the context of quantum diffusion of particles in dynamical disordered 
continua. In this work we will apply these methods to our model of space-time fluctuations and show that 
also in our case one obtains modified wavepacket dynamics dependent on the properties of the stochastic 
model. 



2 Modified quantum dynamics 

In our model space-time is regarded in the low energy limit as a fluctuating entity which consists of a 
classical fixed background on which Planck scale fluctuations are imposed. We assume that they appear 
as classical fluctuations of space-time and model this as perturbations of the metric up to second order 

g^^(x,t) = + hf,^{x,t), (1) 
5^'^(x,<) = ,r -h'^''{^,t) + h^-'i^,t), (2) 
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where <C 1. In this paper greek indices run from to 3 and latin indices from 1 to 3. The second 

order perturbations are given by h^'^ = rj^xh^'^h^^ where indices are raised and lowered with -q^^ and rj^^'^ . 
We consider a scalar field 4>{x) non-minimally coupled to gravity as given by the action 



where g — — det^^j^ is the determinant of the metric and x = (x, t). Here R{x) is the Ricci scalar and 
^ a numerical factor. Variation of the action with respect to (j){x) yields the Klein-Gordon equation 
non-minimally coupled to the gravitational field 



g^-^D^d^^x) - (^-^ + ^Rix)j c^{x) = 0, (4) 

where is the covariant derivative based on the metric g^i/. 

There are three values of ^ which are of particular interest [17]: (i) ^ = is the minimally coupled 
case, (ii) ^ = ^ is required by conformal coupling and can also be derived from the requirement that the 
equivalence principle is valid for the propagation of scalar waves in a curved space-time [18!. (iii) S, = \ 
originates from the squaring the Dirac equation. 



2.1 Non— relativistic limit 

In the next step we calculate the non-relativistic approximation of this equation by expanding the wave- 
function (/)(x, t) = g.i'S{yi,t)/h powers of c? 

5(x, t) = 5*0 (x, t)c^ + 5i(x, t) + 52 (x, i)c-2 + . . . , (5) 

according to the scheme worked out by Kiefer and Singh [1^] • To orders and in the series expansion 
of the Klein-Gordon equation the results are identical to the case where minimal coupling is assumed 
[5] . For obtaining the Schrodinger equation without relativistic corrections in the case of non-minimally 
coupling we give a sketch of the calculations which are done here analogously to [3] . First one has to derive 
the equation of motion to order c°, where non-hermitian terms appear in the Hamiltonian. In order to get 
rid of them one can choose a modified scalar product in curved space or alternatively keep the usual 'flat' 
euclidean scalar product and transform the operators and the wavefunction accordingly [20] . We proceed 
with the latter possibility and finally arrive at 

2m V / 2m 2 

+\{ihd..g'°]^l^. (6) 

Here Acov represents the Laplace-Beltrami operator, {•, •} is the anticommutator, '^^'5 is the determinant 
of the 3-metric gij and g^^ includes the Newtonian potential U (x) . The Hamiltonian is manifest hcrmitian 
with respect to the chosen 'flat' scalar product. 



2.2 The effective Schrodinger equation 

In [5] it is assumed that the particle described by the modified Schrodinger equation has its own finite 
spatial resolution scale. Consequently, only an averaged infiuence of space-time fluctuations can be de- 
tected. This is quantified by means of the spatial average of the modified Hamiltonian over the particle 
scale leading to an effective Hamiltonian 

^ - - {{S"' + o.''{t)) d.d,) - mC/(x) , (7) 

where the tensorial function a*^ (i) consists of the spatial average of squares of the metrical fiuctuations. 
The tensorial function a*-' {t) can be split into a time-average part and a fluctuating part 

a'\t) = + -i'^{t) with (7*^'(t)) = 0. (8) 

It has been schown that the part a*^ leads to a renormalized inertial mass implying an apparent breakdown 
of the Weak Equivalence Principle [3] whereas the fluctuating part leads to a decay of coherences in the 
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energy representation [7j. If we include non-minimal coupling then a spatially averaged Ricci scalar {R)y 
appears leading to an additional term in the perturbation hamiltonian 

Hp = -^{Y'd.d,+^{R)y), (9) 

which commutes with the density operator p of the effective master equation in [7j. Consequently, this 
extra term does not change the decoherence properties of quantum systems. 

In the following sections we will derive the modifications of the dynamics of wave packets as induced 
by the fluctuating space-time metric. This requires an appropriate approximation of the kinetic term. 

2.3 Approximation of the kinetic term 

Now we consider the full hamiltonian ^ and analyze the kinetic term which is composed of the Laplace- 
Beltrami operator and the fluctuating metric quantities. It leads to the expression 

((3)5V4)A_((3)g-l/4) = gV (^((3)gl/4)a^a^.((3)g-l/4)_ Ig^l^y^a^.l^y^' 

-^d^g'^dj In ^((3)^ + g'^d.d, + d.g'^d,. (10) 

All contributions of the kinetic part as well as those terms in the anticommutator of equation ([6]) which 
do not include derivatives acting on the wavefunction give a modified stochastic scalar interaction term 



nx,t) ^ -|-(ei? + g^^(((3)5i/4)a,a,((3)g-V4)) 

(^-\d. In ^^)d, In - \d.g^=d, In 

+ 'p.9'° (11) 
appearing in the Schrodingcr equation 

ihdti> = -^{g'HK,t)d,d,+d,g'^dj)i; + ihg''>d,^ + Vi^,t)^, (12) 

where for simplicity the Newtonian potential g*^" is not included henceforth. Our conclusions will not be 
affected by taking a Newtonian gravitational field in to account. 

Now we will employ some approximations which allows to use established methods for the calculation 
of wavepacket dynamics. We assume that space-time fluctuations act on a space-time scale which is 
much smaller than typical length- and timescales on which the wavefunction ip{:ic,t) varies. Therefore the 
amplitude of the terms containing derivatives of the fluctuating quantities (/i'-' and the trace h) dominate 
in the Schrodinger equation ()12p compared with the magnitude of the other fluctuating quantities. This 
allows us to state that the main contributions from spacetime noise stem from terms containing products 
of first derivatives of the fluctuation quantities . In short notation 

\dhdh\ > {h,h)dd, \dhdh\ > hdhd and \dhdh\>dhd, (13) 

where for simplicity we omitted indices. Therefore we can neglect the terms 

In addition, we note that in the following derivation of the averaged master equation only second 
moments of V will appear, see Eq. ()22p . This would lead to third and fourth order terms if we incorporate 
second order fluctuation terms (third order terms vanish in the average because of the gaussian property 
which we will specify later). Therefore we can neglect second order terms already on the level of the 
Schrodinger equation leading to 

ihdt^l^^~ — l^i}j + V{ic,t)^, (14) 
2m 

where A is the Laplace operator in cartesian coordinates and where the random scalar interaction (jlip 
now reads 

VM = -^(^^R + 5^j((i^)g^/^)d,d,e)g-^/^))). (15) 
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To first order in the perturbation metric this reduces to 

vi.,t) . -^X^R-l^h), (16) 

where h is the first-order spatial part of the trace tr(gpy). To first order the Ricci scalar is given by 

R{x,t) = didkh"'{x,t)- Ah{x,t), (17) 

Note that only spatial components are present which is an effect of the non-relativistic approximation ([5]) . 
Compared to the second term in the fluctuating potential (|16p the non-minimal coupling term has the 
same structure and the same order of magnitude. Thus this contribution to the Schrodinger Hamiltonian 
and to the analysis of space-time fluctuations cannot be neglected in the context of our model and must 
be taken into consideration in the following calculations. 



2.4 Stochastic properties 

Now we have to specify the statistical properties of the fluctuating potential ^^(x, t) as we like to derive 
an effective equation. Owing to the fluctuating metric this term will be interpreted as a gaussian random 
function V{:x.,t) dependent on the quantities h^^{yi^t). We will now continue the calculations for the one 
dimensional case as in the next section the effective master equation is derived for this specific case and 
from now on the spatial coordinate is denoted by x. This yields for the fluctuating potential 

V{x,t) = (^ea>ii(:r,i)-a2^n(x,<) J)) (18) 

= ^dlh,,{x,t). (19) 

As the quantities of interest are the metric perturbation terms, we choose the following statistical condi- 
tions 

{h^i{x,t)) = (20) 
C{x,x',t,t') = {hii{x,t)hii{x',t')) = hl5{t^t')C{x-x'), (21) 

where /iq is the strength of the fluctuations of h{x,t). It has been shown by Heinrichs [TH [T5] that the 
dynamics of wavepackets are unaffected by temporal correlations of the random function h{x,t) if they 
are sufficiently small. This also applies to our model of space-time fluctuations because it is commonly 
believed that quantum gravitational induced fluctuations appear on typical scales given by the Planck time 
Tp and Planck length Ip. Therefore we have chosen a J-correlation with respect to time. However, at the 
moment the spatial correlation function C(x — x') will be left unspecified. For the following calculations 
it is convenient to introduce the correlation function 

{V{^,t)V{^',t')) - Vo'S{t-t')g{x--x'), (22) 

where g{x — x') = 9^9^, C(a; — x') and we defined Vq = g^^o- 

As this is a model for space-time fluctuations and because of the lack of a complete understanding of 
the microscopic structure of space-time one has the freedom to characterize the statistical properties of 
the terms h with a variety of stochastic models given here by the correlator C{x — x'). 



3 Dynamics of the density matrix 
3.1 Effective master equation 

Most physical quantities of interest - in our case the mean square displacement - can be conveniently 
calculated in terms of the reduced density matrix {p{x',x,t)) and its evolution equation. For the sake of 
simplicity we restrict all the calculations to the one-dimensional case. In the following we pursue closely 
the strategy and calculation techniques realized in [16|, I21| , which lead to an exact solution of the master 
equation and thus for the mean-squared displacement of a wavepacket. In addition to this we clarify some 
derivations and assumptions made in _16J and apply the results to our model of space-time fluctuations. 
We define the density operator 

p{x',x,t) =ijj*{x',t)ip{x,t) (23) 
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and set up the master equation according to 

dtpix',x,t) = dti^*{x',t))ijj{x,t) +ij*{x',t)dMx,t) (24) 

= ^ - ^) '^(^'' ^' + ^ (^(^'' ^) ~ ^(^' *)) '°(^'' ^' (25) 

The formal solution of this master equation reads 

pix', X, t) = p{x', < = 0) + ^ dt' (^-^ - ^) Pix',x, t') 



+- dt' {V{x',t')^V{x,t'))p{x',x,t'). (26) 
" Jo 

In the next step we calculate the effective equation of motion by applying the average over the fluctuations. 
Special care must be taken for the terms which are products of the fluctuating potential V{x,t) and the 
density operator. This is due to the fact that the average of the product {V{x, t)p{x' ,x, t)) generally does 
not factorize because the density operator p{x' ,x,t) is a functional of the fluctuation potential V{x,t)^ 
hence p\V]. 

We can handle this by using the fact that we have chosen gaussian fluctuations (see the statements 
preceding Eqns. ^ and A gaussian process is characterized in functional form by the Novikov 

theorem [2? valid for processes obeying {V{x,t)) = (here: {hii{x,t)) = 0), which yields in our case 

{V{x,t)p{x',x,t)) = j dt" j dx"6{t~t")g{x-x")(^-0^^, (27) 
{V{x',t)p{x\x^t)) = Jdt" Jdx"Sit-t")gix'-x")(^^^^y (28) 



where g is the correlation function introduced in equation 

The Novikov theorem follows from a functional Taylor series expansion of p[V] and utilizing the gaussian 
property of hii{x,t) which leads to a factorization of moments j23|. The functional derivative can be 
calculated using the formal solution of the master equation ((26|) and yields 

im^) ^ \[dt'p{x',x,t'){5{x-x")6{t'-t")-5{x'-x")6{t'^t")) 

= ^p{x\x,t"){5{x-x")-5{x' -x")). (29) 



This result used in (|27p yields 



{V{x, t)p{x', X, t)) = ^ y dx"gix - x") {p{x\ X, t)) {5{x - x") - 5{x' - x")) 

= '-^{p{x',x,t)){g{0)-g{x-x')), (30) 

and analogously, from equation 

{V{x',t)p{x',x,t)) = '^{p(:,'^:,,t)){gix-x')-giO)). (31) 

Combining both results and inserting this into the averaged Equation ([25)1 we obtain the effective 
master equation for the reduced density matrix {p{x' , x, t)) 

-(p(.,.,i)) = -_(^___j(p(.,.,t)) 

-^{g{0)^g{x~x')){pix\x,t)). (32) 
3.2 Solution of the effective master equation 

In order to get rid of the time dependence we perform a Laplace transformation according to 

r{x',x,s) = C[{p{x',x,t))]^ dt{p{x' ,x,t))e-'\ (33) 



5 



where the Laplace variable s satisfies Re(s) > 0. The time derivative becomes 

C[dt {p{x\ X, m = sC[{p{x\ X, t))] - {p{x', X, t = 0)) (34) 
and the Laplace transformed master equation reads 

' —^]r{x',x,s)+[s + -^ig{0)-gix-x'))]r{x',x,s) = {pix',x,t = 0)), (35) 



2m \ dx'^ dx^ J ' ' \ h? 

where {p{x' , x,t — 0)) is subject to initial conditions. For practical reasons we introduce new coordinates 
according to X = x + x' and Y = x — x' which yields 

-AX,Y,s)+(s + %{g{Q)-g{Y))\r{X,Y,s)^{p{X,Y,t^Q)). (36) 



m dXdY ' ' ' ' \ 

This is the partial differential equation (PDE) which has to be solved. 

By using the averaging techniques to derive the effective equation (j36p we converted the stochastic 
PDE ((25|) to an 'ordinary' PDE. Note that the information about the fluctuation scenario is encoded in 
the correlation function g{Y) and the amplitude Vq. 

We rewrite the master equation p6p as an ordinary, inhomogeneous first-order differential equation by 
applying a Fourier transformation w.r.t. the variable X 

/oo 
dXe-'^^r{X,Y,s). (37) 
-OO 

This gives 



dvRiK, Y, s)-^(^s+^ (g(0) - g(Y))^ R{K, Y, s) = -^R{K, Y,t^O), (38) 



where R{K,Y,t = 0) = {Tx (p)) {K,Y,t = 0) is the Fourier transform of the inhomogeneous part of the 
ordinary differential equation. 

This equation can easily be integrated and yields 

R{K, Y, s) = exp + G(y)) RiK, Yo, s) - exp + G(r)) x (39) 

m , , f msY'\ , , , ( mG{Y'\ 

^ j^zi \ r,i i 



I ^>^'exp(-^)wr',t = 0)exp(- 



hK ]q ' \ 2hK J V 2fijs: 

where the term R(K, Yq, s) = R{K, Y = 0,s) represents the initial value of the solution R{K, Y, s) and we 
defined G{Y) = f dY' (g(0) - g{Y')). 

3.3 Mean squared displacement 

The mean squared displacement of the particle is given by 



crlit)^ / dxx^ lim {p{x',x,t)) (40) 
and can be expressed in Fourier space as 

Therefore the quantity of interest which has to be calculated is the initial value of i? at F = 0, corre- 
sponding to a; = a;' and hence to the diagonal element of the density matrix p{x',x,t). 

The solution R{K, Y, s) is a quadratically integrable function and must vanish for Y oo. Therefore 
the r.h.s. of equation (|39p must be zero in this limit leading to 



R{K,Yo,s)^— dY'exp(^^^\R{K,Y\t^O)exp 



mGiY'] 
2hK 



(42) 
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We make the substitution r = 2^1?^ (the dimension of this quantity is in fact that of time) and note that 
the r.h.s. represents a Laplace transformation 



/•oo 

R{K,Yo,s) ^2 dTCxp{-sT)R{K, 
Jo 



2hK 



•,i = 0)exp(-G(r)), 



(43) 



if T is positive giving K —t \K\. This equation can be inverted leading to the solution 

^2 rr r ^2h\K\ 



i?(if,ro.r) = 2i?(i^,^fflT,< = 0)expf— ^ f dr' 



5(0) 



(44) 



Inserting this expression into equation (j4ip we get for the mean squared displacement 

cjl{t) = ie-«(^)i?o(y) {dlG{Y) ~ {dkG{Y)f + 2dkG{Y)dk\nRo{Y) - dllnRo{Y) - {dk\nRo{Y)y 

where Y = 2^r and Ro(Y) = R{K, '-^t, t = 0). 



K=0 

(45) 



Note that the unperturbed solution is given by the last two logarithmic terms, where for K — we get 
exp {—G{Y))Rq{Y) — 1, if the initial state i?o is given by a gaussian wavepacket. The information about 
the fluctuation scenario is given by the correlation function appearing in G(Y). In our case we regard 
even correlation fimctions G{Y) = G{—Y) excluding a preferred direction in space. For our purposes it is 
assumed to have the form G(l"2"), for n € N avoiding the introduction of the modulus of the argument 
Y as the correlation function must be differentiable and analytic at the point Y = 0. In this case the first 
derivative dkGiY) vanishes in the limit K — therefore modifications to the free (unperturbed) mean 
squared displacement come from the second derivative df.G{Y). 



3.4 Gaussian correlator 

With this result (j45p we can now calculate the mean squared displacement of a wavepacket. In the 
preceding discussion we imposed several restrictions and assumptions for the argument of the correlation 
function CiY). For this reason and as a first starting point we choose a gaussian correlator 

C(r) = -^expf-^Y (46) 



where a represents a finite correlation length and we choose for the initial conditions a gaussian wavepacket 



R{K, r, t = 0) = Tx 



1 f X^ + Y' 



(27r)i/2a ^'^'^ V 8a2 
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(47) 



= 2 exp ( - ( — + 2K^a^ ) ) . (48) 



We set Y — 2-^^^^t and the mean squared displacement yields 

-'(^) = ^"M + T^Juv/ + -T^^'"- (49) 

4TO^a-^(0) VSTrm^a" 

The first two terms correspond to the expression that is obtained from the free Schrodinger equation 
whereas the last term accounts for superdiffusive behavior. The asymptotic behavior for i oo is 
dominated by the dependence 

al{t) ^ alio) + -=^t\ fort^oo. (50) 



We can state that for long evolution times the wavepacket undergoes superdiffusion if the spatial correlation 
properties of space-time fluctuations are given by (PH)) . The effect of this special model of spacetime 
fluctuations is such that the evolution of quantum particles is superdiffusive. 

The inclusion of a white noise scenario by setting C{x — x') = S{x — x') is not feasible in this framework. 
The reason is that our equations are only valid for analytical correlation functions (analytic at i<' = 0) but 



^The nomenclature is at follows: For the mean-squared displacement {x^) oc the exponent u = 1 renders diffusive 
motion, = 2 is a ballistic motion and i/ = 3 denotes superdiffusion. 
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singular quantities like S{Q) would appear in the expressions which were derived here. These expressions 
would render infinite momentum contributions which must be limited by a momentum-cutofi^. In our 
framework this is motivated by a fundamental minimal length-scale which is given by the Planck length, 
since we consider possible effects emerging from quantum gravity. This leads to the desired momentum- 
cutoff and hence the Dirac ^-function is not the appropriate choice in this context. 

However, a modified behavior for wavepackets was already sketched in [7] where a white noise scenario 
for the fluctuating phase of the wavefunctions was considered yielding 

alit) = alio) + + + 4^P*- (51) 

TO TO^ TO^ 

The quantities appearing here are given by apx{t) = {px + xp)t — 2{p)t{x)t and ap{t) = (p^)f — (p)^ which 
is constant in time. The constant Tp appearing in the dissipative term is the Planck time and given by the 
amplitude of the white noise scenario whereas the average is deflned by the trace over the effective density 
operator {A)t = tr{A {p{t))}. However, for large times t ^ Tp the quadratically increasing ballistic term 
dominates and the dissipative linear term can be neglected thus reproducing the behavior of a wavepacket 
described by a free Schrodinger equation. 

Therefore the result ()50p derived here by introducing a finite correlation length (gaussian correlator) 
represents a new effect emerging from our model for spacetime fluctuations. 



3.5 Other correlators and higher order moments 

In general it is not possible to distinguish different fluctuation scenarios on the level of wavepacket dynamics 
only by the second moment (|¥T|) . This does not suffice to fully characterize the modified wavepacket 
evolution what we demonstrate here by selecting the correlation function 

C(r)=C4Cxpf-^Y (52) 



4a4 J 

where C4 is a normalization factor. We arrive at 

'^^(^)='^^(0) + (53) 

showing no superdiffusion and no deviation from the usual wavepacket dynamics. 
Considering the case C(Y) = Ce exp leads us to 

'^lit) = alio) + -^YTTt/ + ^IT^t'^ (54) 
4to'^(t|(0) v27rTO2a' 

where P is the Gamma function. This expression has the same structure and asymptotic behavior like 
equation ((i^ . 

Generalizing these results to correlation functions C{Y) = C„ exp ^— (^-^^^ ^ one yields again 
the expression (l53|) for all n e N>i\ {3}. 

In order to deal with this fact one must consider higher-order moments ^2;^"(i)^ with n G N and n > 1. 

For correlation functions having exponents of order 2n one must include moments up to the order 2n 
to fully characterize the response of the wavepacket to space-time fluctuations, where the higher order 
moments are given by 

{^'''it))^^^R{K,Yo,T)\K=o. (55) 
For example we get for the correlator ([52|) the asymptotic {t 00) fourth moment 

{^\t))^96at{0) + ^^t'. (56) 

Thus, in order to be able to distinguish between different stochastic scenarios for spacetime fluctuations 
one must characterize the wavepacket spreading by the inclusion of higher order moments. 
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3.6 Connection to holographic noise 

As a side remark we want to argue whether and in which way holographic noise would affect the wavepacket 
propagation. While simple estimates lead to a linear scaling of the measurement uncertainty with the 
Planck length Ip, effects from space-time fluctuations in a holographic fluctuation scenario are amplified. 

The holographic principle states that the number of degrees of freedom of a region of space-time is 
bounded by the area of the region in Planck units Ip. Alternatively we can state that although the world 
appears to have three spatial dimensions, its contents can be encoded on a two-dimensional surface, like a 
hologram. According to this principle length measurements of a distance I exhibit an intrinsic uncertainty 
given by 61 > {llp)^^^, where the exponent 2/3 for the Planck length Ip is characteristic for holographic 
noise H. It can be shown [2S] , [M] that this translates to fluctuations of the metric 

6gt.. > {Ip/lf'^a^,, (57) 

where a^^ is a tensor of order 0(1). This corresponds to a power spectral density (PSD) of 

S{f) = /-^/^(4)V3, (58) 

where / is the frequency and c is the speed of light. We write this in terms of a PSD in momentum space 
obtaining 

S{k) = fc-5/6;2/3^ (59) 

where k is the wavenumber. It can be readily checked that this expression leads to 61 = (^^p)^^'^ by 
calculating {61^ = /i'/7°"^ dk {S{k))^ , where kmax. is subject to a cut-ofF given by the Planck length Ip and 
the length I is given by the experiment. This shows that - by means of a correlation function - a finite, 

2/3 

nonvanishing correlation length is involved where the correlator scales with Ip modifying the wavepacket 

2/3 

evolution accordingly. The Ip -dependence distinguishes the holographic noise scenario from the other 
fluctuations scenarios discussed in the previous sections. 

3.7 Experimental suggestions 

For the experimental observation of the predicted effects on wavepacket dynamics it is necessary that 
long evolution times of quantum particles can be achieved. So far, it seems that no dedicated experiment 
concerning the spreading of quantum mechanical wavepackets has been performed. We suggest to test this 
effect. In this context it is feasible that the gravitational field of earth is suppressed as best as possible 
to facilitate the free evolution of wavepackets. This goal can be reached by performing experiments in 
microgravity environments which is realized for example by the QUANTUS [Mj and PRIMUS projects. 
The first already demonstrated that a BEC (Bose-Einstein-Condensate) in free fall can be realized with 
propagation times of about one second. In the PRIMUS project BECs will be used to realize a matter 
wave interferometer under microgravity conditions. We suggest that in both projects cold atom gases 
could in principle be used to study the long time behavior of wavepackets with high sensitivity. 

4 Summary and conclusion 

We have calculated the mean squared displacement for a quantum particle subject to space-time fluc- 
tuations. As we have seen the properties of wavepacket dynamics are strongly dependent on the spatial 
correlation function. For two special cases we encounter superdiffusion of the wavepacket for the second 
moment on large timescales (characterized by the exponent of the time variable). For other correlators 
the second moment is unaffected and the unperturbed expression is recovered. Therefore higher-order mo- 
ments must be included in the analysis showing modifications of the unperturbed case. This leads also to a 
superdiffusive behavior and - in addition - allows to distinguish different stochastic models experimentally. 
For a white noise scenario we encounter a diffusive behavior for short timescales. On large timescales, 
compared to the typical evolution time of a wavepacket, the unperturbed case is recovered as the ballistic 
term dominates. Furthermore we discussed the influence of holographic noise on the wavepacket evolution 
and argued that it should be dependent on the Planck length appearing with the exponent 2/3. 

Therefore, interesting effects representing deviations from the usual wavepacket dynamics can only be 
obtained by inclusion of general correlation functions having finite correlation lengths. In our approach 
the modifications on wavepacket dynamics depend only on the statistical fluctuation properties and the 

^Here the Planck length is assumed to be the fundamental minimal length. 
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magnitude, whereas the latter is determined by the fluctuation strength independently of the statistical 
model. 

Furthermore, we observe that in our model the strength of the modifications to the free Schrodinger 
equation for the minimally and the non-minimally coupled scenario have the same order of magnitude. 
However, the results for the one-dimensional case are independent of the non-minimal coupling term. For 
the 3-dimensional case we expect modifications coming from the non-minimally coupling terms but they 
will not affect our conclusions significantly. Since the structure of the equations remains the same and the 
form of the correlation functions are identical to those in one dimension we should not encounter major 
differences for the expressions for the wavepackat dynamics. The magnitude and the temporal behavior 
of the modifications should remain the same. Nevertheless, possible effects from anisotropic space-time 
fluctuations cannot be described in the one-dimensional case, of course, and therefore our conclusions do 
not include this possibility. 

We conclude that a possible experimental verification of our model of spacetime foam is feasible by 
means of experiments with cold atoms with long propagation times - preferably in a microgravity envi- 
ronment. This gives the opportunity of setting bounds for different scenarios of space-time fluctuations. 
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